In search of the electron electric dipole moment: relativistic correlation calculations of 
the P,T-violation effect in the ground state of HI"*" 
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We report the first results of ab initio relativistic correlation calculation of the effective electric field 
on the electron, EeB, in the ground state of the HI"*" cation. This value is required for interpretation 
of the suggested experiment on search for the electron electric dipole moment. The generalized 
relativistic effective core potential, Fock-space relativistic coupled cluster with single and double 
cluster amplitudes and spin-orbit direct configuration interaction methods are used, followed by 
nonvariational one-center restoration of the four-component wavefunction in the iodine core. The 
calculated value of EcS by the coupled cluster method is EcS ~ 0.345 x 10^'*Hz/e-cm. Configuration 
interaction study gives E^s = 0.336 x 10^*Hz/e-cm (our final value). The structure of chemical 
bonding and contributions to EcS in HI"^ is clarified and significant deviation of our value from that 
obtained in Ravaine et al. Phys. Rev. Lett. 94, 013001 (2005) is explained. 



Introduction. It is known P, Q that existence of the 
permanent electric dipole moments (EDM) of the ele- 
mentary particles violate two fundamental symmetries: 
space parity (P) and time reversal (T). Considerable ex- 
perimental efforts invested recently to the search for the 
electron EDM de (see 0; S IE IS 0) are primarily con- 
nected with the high sensitivity of d,, to the "new physics" 
beyond the Standard Model (see 2] and references). Po- 
lar heavy-atom diatomics with nonzero projection of the 
total electronic momentum on the molecular axis (see 
below) are currently considered as the most prospective 
objects to search for de because of the great value of the 
effective electric field acting on the unpaired electrons 
in the ground or excited states of such molecules H, . 
The only nonzero component of the effective electric field 
in polar diatomics is the one directed along the molec- 
ular axis and is traditionally written as Ecs = Wd\^\ 
where $7 is the projection of total electron angular mo- 
mentum J to the molecular axis, (see Eq. © and [illi3 
for details). Calculation of the value is needed for 
any experiment on the d^ search using molecules. The 
theoretical methods developed recently (see ^| and 
references therein) allow one to calculate E^s for any di- 
atomic molecule of practical interest with required accu- 
racy, even for such complicated system as excited states 
in PbO [Illli. 

Recently, the EDM experiment of new type, on the 
molecular HI+ cation in a trap, was suggested by Stutz 
and Cornell ^E|- The estimates for the value of E^s 
in two molecular cations, HBr+ and HI+, were made by 
Ravaine et al. [l^ and markable influence of the chemical 
bond nature was emphasized. In the present article, we 
report the first results of ah initio calculation of £^eff for 
the ground state of HI+, leaving out discussion about 
sensitivity of the suggested experiment. 

Following ^3 , Ravaine et al. modelled electronic struc- 
ture of HI+ in [i3 by two limiting approximations: 
"ionic" and "covalent" , where "ionic" approximation cor- 



responds to a neutral iodine atom electrostatically per- 
turbed by a proton. The latter is located at the experi- 
mental equilibrium distance from the iodine nucleus de- 
termined for the HI"*" ground state. The "covalent" limit 
corresponds to the 1+ ion perturbed by the electrostatic 
field from the dipole of the polarized neutral H. We are 
using inverted commas to mark the approximations since 
from the traditional point of view they both correspond 
rather to a van der Waals-type interaction between ion 
and induced dipole. It was found in jl6l | that the value 
of i?off is greatly changed (about six times) depending on 
the approximation made, either "ionic" or "covalent" . 

Model consideration. We calculated molecular dipole 
moment D of the ground HI"*" state with the molecular 
axes origin at the iodine nucleus using the restricted ac- 
tive space self-consistent field (RASSCF) |il| method. 
The details of that calculation can be found on Il9|. Our 
calculation shows that the highest doubly occupied cr- 
orbital is bonding and most "mixed" one among occupied 
orbitals. It is formed mainly by the iodine 5po and hydro- 
gen Is orbitals, where subscript denotes the projection of 
the angular momentum on the molecular axis. Though 
this is not the highest by energy from the occupied or- 
bitals, it gives 77% of the calculated electronic part of 
£), Del — —2.610 (the contribution from the hydrogen 
nucleus is, obviously, l-i?e = 3.08; here and below we 
use atomic units unless the opposite is stated explicitly). 
On the other hand, the valence 7r_i and tt+i orbitals are 
formed mainly by the iodine 5p_i and 5p+i orbitals (mix- 
ture of the 5pi/2,±i/2, 5p3/2,±i/2, and 5p3/2,±3/2 spinors) 
and their contribution to D^i is about 7%. 

The distinctions of the "covalent" and "ionic" approx- 
imations from our consideration can be illustrated on 
the onc-configurational model of the chemical bond in 
III+. The leading (SCF) configuration of this molecule 
in the ground ^113/2 state with Q, — +3/2 (having weight 
0.9 in the correlated wavefunction) can be presented 
as [. . . ]cr^(7r^i_Q,7r+i,^)7r;^]^^^, where a and (3 correspond 
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to the spin projections +1/2 and —1/2, the degenerate 
TT-i^a and 7r+i^/3 spin-orbitals constitute a sheU in the 
used relativistic classification, and the unpaired tt'^i ^ 
state is mainly the 5p3/2,+3/2 spinor of iodine (we distin- 
guish it' from TT just to emphasize their different behavior 
at the iodine nucleus). Let us consider contributions of 
the bonding spin-orbitals Iffa,/}) = C/ki,/?) + dnlcr^p) 
to the hyperfine structure (HFS) properties at the iodine 
nucleus. The occupied shell can be presented as 

+ CiCniaiaf, +a^a'p) (1) 

where the operator of asymmetrization is omitted, \(t^) 
is mainly 5po-orbital of iodine, \cr^) is mainly Is-orbital 
of hydrogen, and Cf -I- C|f — 1 (assume for simplicity 
that |(T^) and \a^) are orthogonal and Ci,Ch are real; 
C/«0.83, C//«0.56 in our calculations). In the spin-orbit 
representation, the atomic 5po orbital can be approxi- 
mately presented as a combination of the 5pi/2 and ?>p^/2 
spinors of iodine with weights 1/3 and 2/3, respectively. 

The contribution in the first line of Eq. |^ corresponds 
to the "ionic" model in [11 [13 (when \Ci\ - 1), two 
terms in the second line can be compared to their "co- 
valent" model (the maximum |C/C/f|=l/2 is, obviously, 
attained for |C/| = |C//| = l/v^) and the term in the 
third line corresponds to the conventional ionic model 
I"^"*"-!!" {\Ch\ ~ 1) that is not considered there. Both 
terms from the second line equally contribute to the elec- 
tric quadrupole HFS constant on iodine (dependent on 
the space-inhomogeneous part of electronic density with 
respect to the iodine nucleus) whereas their contributions 
are completely compensated for the spin-dependent mag- 
netic dipole HFS constants as well as for E^s (see be- 
low). The latter is a consequence of the fact that closed 
shells do not influence on those spin-dependent proper- 
ties; their contributions can become nonzero only when 
polarization and correlation effects are considered. So, 
the only open shell ir'j^^ ^ (5^3/2,-1-3/2) should be consid- 
ered when calculating spin-dependent properties within 
our simple one-configuration model. The weight of the 
rest configurations (obviously, orthogonal to the leading 
one and accounting for correlation without any restric- 
tion on the occupancy of TT±i^cn '^±i,i3, f^'^j ^^p^ etc.) is 
only 0.1 (as is obtained in the calculations discussed be- 
low). Even if one suggests that all the correlating config- 
urations contain singly-occupied Is-orbital of hydrogen 
and a" is also Is, the maximal weight of the HI"*" config- 
urations of type [. . . ]ils\j is smaller than 0.7 that can be 
compared to the weight 1.0 in the "covalent" model of 
Chanda et al. ^3 and Ravaine et al. Moreover, 
in the "covalent" model by Ravaine et al. the contri- 
butions with the weights 2/3 and 1/3 (see Eq. (13) in 
[IG]) in the front of the wavefunction terms containing 
(j^ and (T? (Is of hydrogen) are fixed in accord to the 



spin-coupling rules (with the lowest lying ^n2,i states of 
I"^ and proper dissociation limit) and not varied. This in- 
duces a large artificial asymmetry in contributions of the 
5pi/2,±i/2 states when calculating E^s within the "cova- 
lent" model both at the one-configuration and correlation 
levels. Such asymmetry could be attained in the corre- 
lation calculations only if the configuration in which the 
singlet (T^ pair replaced in the leading configuration by 
the triplet cr^crjj state would have weight « 0.1, thus 
leaving nothing to other correlation configurations. 

Effective P,T-odd Hamiltonian. The terms of our in- 
terest for HI+ in the effective spin-rotational Hamiltonian 
may be written following Refs. 0, [i^l . The P,T-odd in- 
teraction of de with the effective electric field is 

Hd = Wd de{J ■ n) , (2) 

where J is the total electronic momentum and n is the 
unit vector along the molecular axis from I to H. In 
slightly different form of Hd is used: 

Hd - Wd{J ■ n) , (3) 

so the value of de appears explicitly in their final result 
for Ecs. The effective operator 

H. = 2de[l^%) (4) 

is used to express the interaction of de with the inner 
molecular electric field E {cr are the Pauli matrices), to 
avoid large numerical cancellation of the terms with op- 
posite signs because of Schiff's theorem [^H^- After 
averaging over the electronic coordinates in the molecular 
wavefunction, one obtains 

Wdn = ^{^n\Y.He{i)\^n) , (5) 

^ i 

where is wavefunction for the X^n3/2 state. 

To check the accuracy of calculating the wavefunction 
in the vicinity of the iodine nucleus we computed the hy- 
perfine constant A\\ (see [l^) and quadrupole coupling 
constant eQqo, where Q = —710 millibarn is quadrupole 
moment of ^^^I [2^, go is electric field gradient along 
molecular axis. Note, however, that the errors in cal- 
culated eQqo and Ecs are not related closely. As 
our recent calculations showed [l3. H^. the error for Ay 
presents rather a lower bound estimate for the E^g error. 
The quadrupole interaction constant is capable to pro- 
vide useful information about space-inhomogeneous part 
of the electronic density near a nucleus. Unfortunately, 
the eQqo value is not a better measure of the calcula- 
tion accuracy of the effective field on the electron than 
the A|| , first of all because it doesn't depend directly (like 
EcS and ) on the electronic spin density near the heavy 
nucleus. In the parameters of the Frosch-Foley effec- 
tive spin-rotational Hamiltonian (3| were obtained for 
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the ground state of HI"*" . The connection of the Frosch- 
Foley parameters a, 6, c to is [ssj : 

A|| = ^(aA+(6 + c)S) , 

where A is the projection of the angular electronic mo- 
mentum on the molecular axis and T,—Sz is z-projection 
of the electronic spin. Accounting for the calculated value 
of G||-factor that is 2.0001 and \n\ = 1.4998, the ground 
state of HI+ can be reliably classified as ^^3/2- 

Methods and calculations. A 25-electron generalized 
relativistic effective core potential (GRECP) ,25] for io- 
dine (its gaussian expansion can be found on our website 
P^J) is used at the first step of the two-step calculations 
of HI+, so that the inner shells of iodine, Is — 3d, are 
absorbed into the GRECP and the 4s, Ap, Ad, 5s, and 5p 
electrons (as well as an electron of hydrogen) are treated 
explicitly. Two calculations are carried out. In the first 
one, only seven external electrons of iodine are corre- 
lated whereas its 4s, Ap, Ad shells are "frozen" within the 
GRECP approach when employing the level-shift tech- 
nique [2^. Thus, a 7-electron GRECP version is, in fact, 
used in the first series of the HI+ calculations. In the 
other calculation, all 25 electrons are explicitly corre- 
lated. The terms with the leading configurations a'^ir^ 
are calculated where a and tt are the highest occupied 
molecular orbitals. The correlation spin-orbit basis sets 
are optimized in atomic two-component relativistic cou- 
pled cluster calculations of iodine with single and double 
cluster amplitudes (RCC-SD) using the scheme suggested 
in H^H^I- As a result, the basis [5s5pid2flg] was gener- 
ated. As our investigation shows removing of g-function 
from the basis set changes the RCC-SD results for 
and Ecff on the level of 1%. Thus, contribution from g- 
function to the calculated values can be negligible and 
the basis reduced to [5s5p3d2/] without loss of accuracy. 
Such iodine basis was used in 25 electron configuration in- 
teraction (CI) calculations of HI+. For hydrjsgen, the re- 
duced [4s3p2c?] correlation-consistent basis [2S| was used. 

The HI''" calculations start from a one-component 
closed shell SCF computation of the ground state of the 
neutral HI molecule using the spin-averaged GRECP for 
iodine. Two-component Fock-space RCC-SD molecular 
calculations or spin-orbit direct CI (SODCI) calculations 
are then performed. 

RCC-SD method: The details on the Fock-space 
RCC-SD method can be found in Ref. 29, 30] and refer- 
ences therein. The program package RCCSD is used in all 
RCC calculations mentioned further in the article. The 
Fock-space RCC calculations start from the ground state 
of HI and use the scheme: 

HI -> HI+ (6) 

with an electron removed from the tt, tt' orbitals. 



SODCI method: Spin-orbit direct CI approach with 
the selected single- and double-excitations from some 
multiconfigurational reference states [U HI] is employed 
on the sets of different AS many-electron spin- and space- 
symmetry adapted basis functions (SAFs). In the SODCI 
code, the double C2V group, C|„, is used to account for 
the spin and space symmetry of the HI"*" molecule, in- 
stead of the more restrictive symmetry group , which 
could, in principle, be employed. In the classification 
scheme, the doubly degenerate ground state has the com- 
ponents only in the irreducible representation E*. 

The SODCI calculations exploiting relativistic scheme 
of configuration selection |3^ start from some space of 
the reference functions: for 25 correlated electrons 4415 
SAFs (see Table IIJ were included in the reference space 
("main" configurations). These SAFs had the largest 
coefficients in the probing CI calculation. The single 
and double excitations from this reference space produce 
about 3 X 10^ of SAFs. Only the most important of them, 
selected by second-order perturbation theory for chosen 
thresholds Ti (see Table ID, were included in the sub- 
sequent CI calculation. About 1.6, 5.7 and 13 millions 
of SAFs were selected for thresholds Ti=0.01, T2=0.001, 
T3=0.0003, correspondingly. 

Since we are interested in the spin-dependent proper- 
ties determined mainly by the electronic wavefunction 
near the iodine nucleus, the shape of the valence and 
outer core four-component molecular spinors are restored 
in the inner core of iodine that is done in the paper 
within the nonvariational one-center restoration scheme 
(NOCR) (see [H El 13113 and references therein). The 
RCC calculation of Ees employs the finite field method 
(see Refs. miilll^). In the SODCI calculations con- 
ventional approach with the density matrix calculation 
for CI wavefunction was used [l^ . 

Results and discussion. The results of the RCC and 
SODCI calculations for 7 and 25 correlated electrons of 
HI+ are presented in Table HI The internuclear distance 
is 3.08 a.u. in accord to the experimental datum n7ll. 

It should be noted that the authors of paper llg con- 
sidered their "covalent" result as the final one and pre- 
sented their "ionic" result only for comparison. The re- 
sults of our RCC and SODCI calculations give essentially 
different EcS value than the one obtained in by the 
configuration interaction calculation for the "covalent" 
approximation. Particularly, the sign of E^g is opposite 
to that by Ravaine et al. One can see that accounting 
for correlations with the iodine core electrons (occupying 
the shells 4s, Ap and Ad) practically doesn't change the 
value of Ees- The importance of accounting for correla- 
tions can be seen by comparing the results of RCC-S and 
RCC-SD calculations. In the RCC-S calculations (only 
with the single-body cluster amplitudes) effect of "spin- 
polarization" is taken into account analogously to "unre- 
stricted" Dirac-Hartree-Fock (DHF) calculation. Inclu- 
sion of electron correlations in the RCC-SD calculation 
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TABLE I: Calculated E^a (in xlO^'' Hz/(e-cm)), yly (in 
MHz) and quadrupole interaction value eQqo (in MHz) for 
the ground state of H^^'^l"'". The one-center expan- 

sion by s,p,d spinors within the iodine core is used in the 
NOCR scheme. Experimental values for is 1021 MHz and 
for quadrupole coupling constant eQqo is —712.6 MHz. 



Method 






^11 


eQqo 


work [16] "ionic" 


approx. DHF 


-0.09 






work fl6] "covalent" approx. CI 


-0.49 






AGRECP/SCF/NOCR calculations 






7 electrons 








resticted SCF 


25 electrons 


0.008 


949 


-647 


resticted SCF 




0.010 


1024 


-667 


GRECP/RCC/NOCR calculations 

7 electrons 




RCC-S 




0.206 


863 


-719 


RCC-SD 


25 electrons 


0.347 


881 


-708 


RCC-S 




0.226 


906 


-807 


RCC-SD 




0.345 


962 


-752 


GRECP/SODCI/NOCR calculations 




Threshold 


SAF number 








(mHartree) 


7 electrons 








Mains only 


7 786 


0.294 


984 


-687 


0.001 


676 397 


0.335 


895 


-711 


0.0001 


1911282 
25 electrons 


0.336 


892 


-709 


Mains only 


4415 


0.333 


1063 


-778 


0.01 


1600 012 


0.299 


975 


-738 


0.001 


5 712 946 


0.329 


971 


-743 


0.0003 


12 678 133 


0.336 


968 


-745 



AGRECP/SCF/NOCR scheme, the value is 0.6 MHz. 

Our results of the SODCI calculation (our final values) 
for the Ecs, and eQqo properties are in close agree- 
ment vifith the RCC-SD values. It means that higher- 
order cluster amplitudes do not contribute largely to 
these properties. Besides, outercore- valence correlations 
practically do not influence on the value of Ecs- On the 
other hand the value of A^^ is increased for about 10% 
when outercore correlations are taken into account. It 
was noticed before that rather good accuracy in the cal- 
culated A|| value gives us just a lower bound for the ac- 
curacy of Ecf[. Taking into account weak dependance of 
EcS from the outercore-valence correlations we estimate 
the accuracy of E^s calculation in 10%. 

In any case our calculations show that the absolute 
value for E^s in X^Il3/2 of III+ is much lower than 
that in YbF, 6.0x10^^ IIz/(e-cm), and in the metastable 
a(l) state of PbO, 6.1x10^^ Hz/(e-cm). Thus, HI+ can 
be perspective candidate for experiments on the EDM 
search provided that the experimental scheme is im- 
proved to reach much better statistics or coherence time, 
than that in on-going experiments on YbF and PbO. 
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changes E^s on about 60%. At the same time value of 
A|| is changed only on 5%, that shows that the structure 
of correlation contributions to A|| and EcS is very differ- 
ent. The same is valid for the eQqo constant, in which 
correlations contribute less than 10%. 

The restricted open shell SCF calculations presented in 
Table m were performed with the spin-averaged GRECP 
(AGRECP) for iodine. The value of E^g is more than or- 
der of magnitude smaller in AGRECP/ SCF/NOCR cal- 
culations than in GRECP/RCC-S/ NOCR ones that in- 
dicates critical importance of accounting for one-electron 
spin-orbit and polarization effects on valence shells in 
calculation of Egg. Similar situation was observed in cal- 
culations on the a(l) state in PbO We would like 
to emphasize that after applying the NOCR procedure 
the proper, four- component shapes of molecular spinors 
in the core of Pb are restored having appropriate rela- 
tivistic behavior at the Pb nucleus both after GRECP 
and AGRECP calculations. In the experiment on III+ 
in the rotating electric field |0| the knowledge of the 
hyperfine coupling value of the proton spin to the molec- 
ular axis can be usefuU. We calculated the value of A|| 
on the H nucleus in the fremework of above-described 
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